DISTRIBUTION OF THE BROWNIAN MOTION ON ITS WAY TO 

HITTING ZERO 
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Abstract. For the one-dimensional Brownian motion B = (Bt)t>o, started at x > 0, and 
the first hitting time r = inf{t > : Bt = 0}, we find the probability density of B U t for a 
u S (0, 1), i.e. of the Brownian motion on its way to hitting zero. 



1. Introduction 

The following problem has been recently addressed in [5], [6]. The authors considered a 
continuous time subcritical branching process Z = (Zt)t>o, starting from the initial population 
of size Zq — x. As is well known, Z t gets extinct at the random time T = inf{t > : Z t = 0}, 
and T < oo with probability one. What can be said about Z T / 2 , i.e. the population size on 
the half-way to its extinction? While the complete characterization of the law of Z u t with 
u = 1/2, or more generally u E (0, 1), does not seem to be trackable, it turns out that under 
quite general conditions 

x u ^Z uT — Cc- U e~ ur] , (1.1) 

x — >oo 

where the convergence is in distribution, C and c are constants, explicitly computable in terms 
of parameters of Z and r\ is a random variable with Gumbcl distribution. 

In this note we study the analogous problem for one-dimensional Brownian motion B = 
(Bt)t>o, started from x > 0. Hereafter we assume that B is defined on the canonical probability 
space (ri,^ 7 , P x ) and let r denote the first time it hits zero, i.e. r = inf{£ > : B t — 0}. 

Theorem 1.1. For x > andu G (0, 1), the distribution P X (B UT < y) is absolutely continuous 
with the density 



Pi^x;y) = — r — ^ 2 TTl _l ^2/ 1 \ _j_ 2 p V- 2 ) 

Tr\(y — xy{\ — u) + y uj{(y + x) z {l — u) + y u) 

Remark 1.2. Notice that p(u,x;y) decays as cx 1/y 2 and hence its mean is infinite. Such 
behavior, of course, stems from the possibility of large excursions of B from the origin, before 
hitting zero. 

Remark 1.3. The formula (|1.2|) implies that x~ 1 B UT has the same law under as B UT under 
Pi, or using different notations, 

^Kr^^Kr^ (1-3) 

where B x stands for the Brownian motion, starting at x > 0, and r(x) = ini{t > : Bf = 0} 
(i.e. -B^/jn has the density p(u, This scale invariance does not seem to be obvious at 
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the outset and should be compared to , where the scaling depends on u and holds only in 
the limit. 

In the following section we shall give an elementary proof of Theorem 11.11 In Section [3] our 
result is discussed in the context of Doob's ^.-transform conditioning. 

2. Proof 

Let 8 > and define^ t$ := 5[t/8\ . Recall that r has the probability density (see e.g. [2]): 
/(x; t) = ^-P X (T <t) = -^=e~ x2/2t , t > 0, x > 0. (2.1) 

Let M s ,t '■= inf s < r<t B r and </>(■) be a continuous bounded function, theifl 

oo 

EM b uts) = ^{B UT5 )l{r G [5k, S(k + 1))] 

k=Q 



= £ x <t>{Bu5k)l(r G [6k, S(k + 1)) 

oo 

= ]T E^(B u5fe )/(M ,ifc > 0,M S k,s(k+i) < 

oo 

J2 E *<i>{Bu6k)l(Mo,uSk > 0)P*(&u6kjk > 0,M 5M(fe +i) < 0|^ 



fe=0 
oo 



k=0 



(^M u Sk,8k > 0,M5fc,5(fe+l) < 0|-^u<5fe 

oo 

)Px(r G [0,5)) + E x </.( J B„5 fe )7^M ,«5fe > o) x 

Pb„« (t G [(1 - u)<Sfc, (1 - u)** + 5)) 



fc=i 



(x) f(x;t)dt + J2E x (j)(B u 5 k )l(Mo,u8k>0) f(B uSk -t)dt 

JO fe=1 v 7 J(l-u)6k 

pS poo f °° /.(5(fc+l) > 

(z) / f{x;t)dt+ / </>(yK V / q(x,uSk,y)f(y;t - u5k)dt\dy 

JQ JQ I k=l^ Sk ' 



pS pOD f pOO \ 

= 4>(x) J f(x;t)dt + J <Kv)\] q(x,u[t/S\6,y)f(y;t-u[t/6\5)dt\dy 

p5 poo ( poo ~\ 

= <f>(x) f(x;t)dt+ <t>{y)\ q(x,u\t/S]S,y)f(y;t + 5-u\t/S]S)dt\dy, (2.2) 



where q(x, t, y) is the probability density of P x (Mo,t > 0, B t G dy) with respect to the Lebcsguc 
measure (see e.g. formula 1.2.8 page 126, [2]): 

q(x,t,y) = (-l=e-^) 2 / 2t - -L e -C»+-) a /»\ ^ > . (2.3) 

1 V 27TI V Z7TI J 



1 La;J stands for the integer part of x G M and := \_x] + 1 
denotes the indicator function 



By continuity of the densities (|2.1I) and (|2.3|) . for any fixed x > and u <E (0, 1), the function 

F S (t, y):=q(x,u \t/S] 5,y)f(y;t + 6 - u \t/S\ 5) 

converges to 

hmF s (t,y)=q(x,ut,y)f(y;t-ut), Vt > 0, y > 0. 

d^O 

In Lemma |2 . 1 1 below we exhibit a function G(t,y), independent of 5, such that 

F S (t,y)<G(t,y), \f{t 7 y)eR 2 + and / G{t lV )dtdy < oo, (2.4) 



-+ 



and hence, the dominated convergence and (|2.2| imply 



lim E x (f>(B UTS ) = Jim / 4>{y)F s (t,y)dydt = \ 4>(y)q(x,ut,y)f(y;t-ut)dtdy. (2.5) 



On the other hand, lim^o t$ = r, P^-a.s. and thus by continuity of B t , lim^o B UT5 = B UT , 
P x -a.s. for any u £ (0, 1). Thus, by arbitrariness of <fi, (|2.5p implies that the distribution of 
B UT has the density: 



p(u,x;y):= / q(x,ut,y)f(y;t-ut)dt. 

A calculation now yields: 

p(u,x;y) = r V - e -v 2 mi-u)i_^ e -(y-x) 2 /2ut _ _^ =e - {y+x f/2ut 1^ 



o V2tF(*(1 -u)f 2 \V2^d v / 2^rf 



2tt(1 - «)3/a«Va y„ *a 
and by a change of variables 



( \ ^ My-z) 2 y 2 X' 1 / (j/ + ^) 2 _jr_V 

P{U,x,y } 27r(1 _ u)3 /2 u i/2|V 2w + 2(1- M )J V 2u + 2(1 -u)J 
2yV^ f 1 1 



27iVl - u \ (y- x) 2 {l -u) + y 2 u (y + x) 2 (l - u) + y 2 u 

\/u(l — u)4xy 2 
7r|(y — x) 2 (l — u) + y 2 it||(y + x) 2 (l — u) + y 2 w| 

The statement of the Theorem 11.11 now follows from: 
Lemma 2.1. (|2.4[) holds with G(t,x) defined in (|2.6p below. 
Proof. Set t 5 := [i/ 5 !^ so that f < ^ < t + 5, and 



< /(f < S) 1 c-fa-^V^-) 2/ e _ly»/(f + (l_ u )J) 



^ p -iy 2 /(t+'5-«*' 5 ) 



'w<5 ((1 - u)<5) 



3/2 



/(f > 5)—!— J e -fe-^) 2 /2«t 5 _ e -(y+ 2; ) 2 /2'«t' S I ?/ e -±y 2 /(t(l-u)+6) 

Vwi\ ) ((t + 8)(1 - u)f /2 

--: I(t < S)A + I(t > 5)B. 
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Since the function z 2 e Cz with C > attains its maximum 4e 2 jC 2 on the interval [0,oo) at 
2 := 2/C, 



y/u(l - u) 3 S 2 { V 2u 2(2-u)y<5 



(y-x) 2 1 y 2 



i/it(l - u)3 V 2u 2 (2 - u) 

Similarly, for t > S, 

B = y p -(y-xf/2ut s L _ e -2xy/«t' 5 ] >e -iy 2 /(t(l-«)+'5) 

v/u(l - u)3i(i + <S)3 I J 

< ^ J_p-(y-^) 2 /2^(t+5) / ] _ e -2xs,/«tL -ij,V(t(l-u)+5) 

" v/«(l - <z) 3 1 2 1 / 

Hence for 5 € (0, 1] we have the bound 

Fs(t, y) < ^=L= ( kzi£ + i-^) " a 7( t < 1) + 
- it) 3 V 2 « 2(2-u)y 



-(y-x) 2 /4ui J 1 _ ^xy/ut I - %y 2 / (t(2-u)) = . QUy\ (2.6) 



^(1 - u) 3 t 2 

Since for u € (0, 1) and x > 0, the quadratic function is lower bounded: 

(y - x) 2 y 2 > ^_ 
u (2-u) ~ 2 ' 

the first function in the right hand side of (|2.6[) is integrable on M. 2 ,. Further, 

y 1 c -(v-xf/iut [ 1 _ e -2 2 ;; / /t l t] >e -iy 2 /(t(2-«)) dt 

o v/u(l - M ) 3 * 2 I J 

(y-x) 2 y 2 V 1 ({y-x) 2 y 2 2xy N 



^/u(l - uf I V 4" 2(2- it) y V 4w 2(2 -u) 
4yit(2 - u) f 1 1 



y/u(l - u) 3 I (2/ - a:) 2 (2 - u) + 2zj 2 u (y - x) 2 {2 - u) + 2y 2 u + 8xy(2 - u) 

32w(2-u) 2 xy 2 /v / u(l-w) 3 
|(y - x) 2 (2 -u) + 2y 2 u}{(y - x) 2 (2 - u) + 2y 2 u + 8xy(2 - it)} 

The latter function decays as oc 1/y 2 as y — > oo and is bounded away from zero, uniformly 
in y > 0, and thus is integrable on R+. Since the last term in the right hand side of (|2.6|) 
is nonnegative, by Fubini theorem it is an integrable function on M 2 ^ for all u € (0, 1) and 
x > 0. □ 

3. A CONNECTION TO DOOB'S ^-TRANSFORM 

In this section we show that the random variable B UT has the same density as the so 
called scaled Brownian excursion at the corresponding time, averaged over its length. The 
latter process is defined by conditioning in the sense of Doob's ^.-transform, and it would be 



natural to identify this formal conditioning with the usual conditional probability. While in 
the analogous discrete time setting, such identification is evident, its precise justification in our 
case remains an open problem. 

For a fixed time T > 0, let R = (Rt)t<T be the 3-dimensional Bessel bridge R = (Rt)t<T, 
starting at i?o = x and ending at zero. Namely, R is the radial par1@ 

R t = \\V t \\, te[o,T], (3.1) 

of the 3-dimcnsional Brownian bridge V — {Vt)t<T with Vo — V and Vt = 0: 

V t = v + W t -^{W T + v), <G[0,T], 

where »gR 3 with \\v\\ = X and W is a standard Brownian motion in IR 3 . 

The law of R coincides with the law of the scaled Brownian excursion process, which is 
defined as "the Brownian motion, started at x > and conditioned to hit zero for the first 
time at time T" . Here the conditioning is understood in the sense of Doob's /i-transform (see 
Ch. IV, §39, [7], and [T], [3] for the in depth treatment). 

On the other hand, one can speak on the regular conditional measure induced on the space 
of Brownian excursions (started from x > 0), given r = inf{t > : B t — 0}. More precisely, 
let £ be a subset of continuous functions C[o i00 )(M), such that for all u> € E, oj(0) = x and 
for each to there is a positive number £(w), called the excursion length, such that u(t) > 
for < t < £(ui) and ui(t) = for all t > £(oj). E together with the smallest er-algcbra £ , 
making all coordinate mappings measurable, is called the excursion space (see §3, [2] for the 
brief reference and p] for more details). Let /j, x (T, ■), T € [0, oo) be a probability kernel on 
the excursion space (E,£), i.e. a family of measures such that T i— > ji x {T,A) is a measurable 
function for all A € £ and (J, X (T, •) is a probability measure on £ for each T > 0. By definition, 
Hx(T, •) is a regular conditional probability of B t /\ T given r, if for any bounded and measurable 
functional F on (E,£): 



E x F(B. At )I(t € A) = F(u>)n x (s,du)f(x;s)ds, VA e 

J A J E 

where f(x;t) is the density of t, defined in (|2.1|) . In particular, for any bounded measurable 
function (f> and some u G (0, 1), 

/•oo /• 

E X 4>(B UT ) = / / 4>(u>(us))n x (s,dLu)f(x;s)ds. (3.2) 
Jo je 

We were not able to trace any general result, from which the identification of ^ x {T,dui) with 
the probability v x (T,dtu), induced on {E,£) by the aforementioned Bessel bridge R, could be 
deduced. While the latter, of course, is intuitively appealing, its precise justification remains 
elusive (some relevant results can be found in [4]). The calculations below show that 

4>{uj(us))v x (s, du>)f(x; s)ds = E x <f>(B UT ), (3.3) 

/o Je 

indicating in favor of such identification. 

For a fixed T > and u G (0, 1), the distribution of R u t, i-e- the restriction of v x (T, duj) to 
the time t := uT, has a density q U T(x;y) with respect to the Lebesgue measure dy, which can 
be computed as follows. We have 

r t(T-t) 



EV t = v(l-t/T), cov(V t )=I- 



T 



denotes the Euclidian norm in '. 
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where / is 3-by-3 identity matrix. Notice that the law of the Bessel bridge R in ()3.1|) doesn't 
depend on the particular v as long as = x and it will be particularly convenient to carry out 
the calculations for the specific choice v = (2, 0,0). Fix a constant u e (0, 1) and let ^1,^2,^3 
be i.i.d. standard Gaussian random variables. Then, for t :~ uT, 



RuT — 



2 



= \j (\/Tu{l - «Ki + x{l - w)) 2 + Tu(l - u)g + Tu(l - u)& 

The random variable 9 := £| + £| has X2 distribution, which is the same as the exponential 
distribution with parameter 1/2 and hence 



d ^ 2 



RuT = bxU + a) +8, (3.4) 



where £ is written for £1 and a := xW (1 — u)/Tu and b := — u)Tu are defined for brevity. 
The density of (£ + a) 2 is given by: 

h{z) := ±P( « + af < z) = ± -±=e-^?l 2 dx = 



dz V / dz y/2~n 

2y/2Hyj~z~ ^ ' ' 

The density of (£ + a) 2 + 9 is the convolution of /1 and the exponential density with parameter 
1/2: 

Mv) ■= [ V fi{z)f 2 (y-z)dz= / y ^^( e -(v^-) 2 /2 + e -(v^ +a ) 2 /2)l -i/2(y-*) dz 
Jo Jo 2V2n^/z \ '2 

1 (z-(*-af/2 + e -(z+af/2y-l/2(y^) d ~ 



2y/2n 

e -a 2 /2-y/2 fjy / _ _ N p -a 2 /2-y/2 fy/y 

2y/2ll 

e -a 2 /2-y/2 



dz = -== / cosh(za)dz 



2na 



■ sinh(^ya). 



Consequently, the density of y/(C + a) 2 + 9 is given by 

h(z) := 2z/ 3 (z 2 ) = V2e _ a ze- z2 / 2 sinh(za), 
V7ra 



and, finally, the density of by/ (£ + a) 2 + 6* is 

1 , / »n V2e-° 2 / 2 _. 



h{z) := -U(z/b) = v 7- , 2 ze--/ 2fc2 sinh(za/6) = e -(-«/6)'/2 _ e -(a+*/»)'/a 

6 y/nab z y/2-Kab 2 



Hence by (|3.4j) . i?„T has the density 

, v V 1 f f 1 (ar(l-u)-J/) 1 

quT(x]y) = j= - —== =i cxp 1 



2nx(l-u) y/Tu{\-u) I V T 2u(l-u 

cxp 



1 (x(i-u) + y y 

T 2u(l-u) 
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We shall abbreviate by writing 

= (x(l -u) -y) 2 = -u) + y) 2 = y 

1 2«(l-u) ' 2 2u(l-u) ' 3 ^FantV2(l - «)3/a 

Then 



<j>(u{ut))v x (t,dw)f{x;t)dt= I <j>{y) I q ut (x;y)f(x;t)dt 

Jo Jo 

DO />OC -1 />00 

m] Cs^{e- Cl/t - e- c ^j^=^e- x2 / 2t dtdy =: J cj ) (y)p(u,x;y)dy 
and by a change of variables, 
p(u, x; y) = C 3 ^= f e -(Oi+x a /a)* _ e^^ 2 / 2 )')*- 2 ^ = 

\/27T ./n ^ J 



C 3 



I Ci + x 2 /2 C 2 + x 2 /2 



A calculation yields, 



(1 - u){x - yf + y 2 u 2 (l-u)(a; + ?/) 2 +r 



Ci + x 2 /2 = ^ a/ ' a - - and C 2 + x 2 /2 



ii 



2w(l-u) ' 2ti(l-u) 

and, consequently, 



p(u,x;y) 



2u{\-u) 2u{\-u) 



2irxu 1 / 2 {l - uf/ 2 v/2tF 1(1- «)(a! - 2^) 2 + ?A (1 - + 2/) 2 + 2/ 2 " 
4a;y 2 y / M(l — u) 



7r{(l — u)(x — y) 2 + y 2 it}{(l — u)(x + y) 2 + y 2 u} 
which in view of (|3^2]) and ([TT2|) . imply ifO]) . 
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